Entanglement in a first order quantum phase transition 



m 
o 
o 

(N 
o 

Q 



■4— » 

o3 



o 
o 



> 
o 
m 

(N 

o 

-I— > 



T3 

c 

o 
o 



X 



Julien Vidal^'Q Remy Mosseri, 1 '^ and Jorge Dukelsky 2, 

'Groupe de Physique des Solides, CNRS UMR 7588, 
Campus Boucicaut, 140 rue de Lourmel, 75015 Paris, France 
2 Instituto de Estructura de la Materia, CSIC, Madrid, Spain 

The phase diagram of spins 1/2 embedded in a magnetic field mutually interacting antiferromag- 
netically is determined. Contrary to the ferromagnetic case where a second order quantum phase 
transition occurs, a first order transition is obtained at zero field. The spectrum is computed for 
a large number of spins and allows one to study the ground state entanglement properties which 
displays a jump of its concurrence at the critical point. 
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One of the most fascinating feature of the quantum 
world is certainly the entanglement which has no clas- 
sical counterpart. Celebrated by the pioneering work of 
Schrodinger yj, and Einstein, Podolsky and Rosen |2| 
about the nonlocality, ubiquitous in the field of quantum 
information 0. Qj. Q entanglement properties of quan- 
tum systems have recently attracted much attention in 
the context of phase transitions. In the various mod- 
els studied such as spin chains in a transverse magnetic 
field Ji Q, II EH Hil, spin ladders [13], spin simplex 
[T3L ITif , Hubbard model 0| the ground state entangle- 
ment has been shown to be strongly modified at the crit- 
ical point raising the question of the universality of these 
behaviours. 

In this Letter, we consider a system where N spins 1/2 
embedded in a magnetic field h that mutually interact. 
We focus here on the antiferromagnetic case, the ferro- 
magnetic one being discussed in Ref. [liij . The symme- 
tries of the Hamiltonian allow us to considerably simplify 
its diagonalization and to determine the phase diagram in 
the thermodynamical limit. A first order quantum phase 
transition is found at zero field whereas in the ferromag- 
netic case, a second order transition occurs at nonvan- 
ishing field. Next, we study the entanglement properties 
of the ground state by computing its concurrence which, 
roughly speaking, measures the two-spin quantum cor- 
relations |l6j . This concurrence which is nontrivial for 
h > 0, is shown to be discontinuous at the transition 
point where it switches to zero. 

We consider the following Hamiltonian which general- 
izes the model introduced in Ref. [T^ : 



i< 3 



= -^( S 2 x +jS 2 y )-2hS z + ±(l +1 ), (2) 

where the cr^'s are the Pauli matrices and S a = ^ c^/2. 
The prefactor 1/N is necessary to get a finite free energy 
per spin in the thermodynamical limit. Without loss of 
generality, we will set h > in the following. The Hamil- 
tonian H preserves the magnitude of the total spin and 



does not couple states having a different parity of the 
number of spin pointing in the magnetic field direction 
(spin- flip symmetry), namely: 



[H,S 2 ] = 0, 



= 0, 



(3) 
(4) 



for all anisotropy parameter 7. In addition, it is straight- 
forward to show that the full spectrum of H is odd under 
the transformation A — > — A and even under h — > —h . 
Furthermore, since H writes in terms of the total spin 
operators, the degeneracy of each eigenvalues belonging 
to a spin S sector is at least equal to the number of spin 
S representations which is simply given by : 



D s = 



N 



\N/2-S 



N 

N/2 - S - 1 



(5) 



for all N. This implies that the full spectrum is ob- 
tained by diagonalizing only one representation of each 
spin sector S which allows us to deal with a large number 
of spins. Denoting by {\S, M)} an eigenbasis of S 2 and 
S z , the matrix elements of H reads : 



(S',M'\H\S, M) 



?s,s> 



A 



N 

2hM 



-N/2 

( s~ S~ £ 
\ a M-l a M °M',M- 

S + S + x 
a M+l a M 0M',M+2 



(1 + 7) {S{S + \)-M 2 
A(l-7) 



8m> 



M'M 



2N 



(6) 



where afj = y/S(S + 1) - M(M ± 1) These expressions 
which are valid for any value of the parameters (A, h, 7) 
and for any N generalize to any spin sector those given 
in Ref. for the Dicke subspace (S = N/2). In the 
isotropic case 7 = 1, one further has [H, S z ] = so that 
H is diagonal in the basis {\S, M)}. 

The antiferromagnetic nature of the coupling between 
spins considered here (A < 0) completely modifies the 
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phase diagram of H as compared to the ferromagnetic 
case. A simple mean- field approach analogous to the one 
presented in Ref. 0, 0] can be performed and predicts 
a first order phase transition in the zero field limit for 
any positive 7. The magnetization (per spin) in the z 
direction of the ground state is given for all 7 > by: 



1(5 2 ) = \sgn(h) 



(7) 



where sgn denotes the sign function which vanishes for 
h = 0. In the thermodynamical limit, the ground state is 
thus the fully polarized state for h 7^ 0. Contrary to the 
ferromagnetic case, the ground state for a finite arbitrary 
N does not, a priori, lie in the symmetric representation 
of the permutation group spanned by the Dicke states 
fl3j |. Indeed, in the zero field limit (h — ► 0), the ground 
state shall belong to the lowest spin sector that minimizes 
the interaction term. The main issue is thus to determine 
the values of h for which level crossings appears. 
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FIG. 1: Ground state energy of each spin sector as a function 
of the magnetic field for 7 = 1/2 (A = — 1). In the even 
case N = 6 (left) all level degenerate at h SUS y (•) whereas in 
the odd one N = 5 (right) a cascade is observed from the 
S = N/2 to S = 1/2 sector. 

For illustration, we have displayed in Fig. ^the ground 
state energy of the different spin sector for a small num- 
ber of spins. Two different scenarios arise according to 
the parity of N. In the even TV case an additional symme- 
try allows us to give a complete description of the ground 
state properties. Indeed, as recently pointed out by R. G. 
Unanyan and M. Fleischhauer j^jj, when 



N 



(8) 



the Hamiltonian is supersymmetric |2JJ . When the SUSY 
condition (jSJ) is fulfilled for an antiferromagnetic interac- 
tion (A < 0), the authors of Ref. [2(| claim that the 
ground state is nondegenerate and given by : 



where A is a normalization constant and where \N/2, y ) 
denotes the eigenstate of S 2 and S y with eigenvalues 
(N/2)(N/2 + 1) and respectively. The parameter rj 
is determined by the relation [2JJ : tanhr; = X/N. 
Even though \ipo) is a ground state of H with total spin 
S = N/2, the spectrum at the SUSY point is however 
highly degenerate. Indeed, at the SUSY point h susy , all 
the lowest eigenvalues of each spin representation are 
equal (Eq = A(l+7)/2) and the degeneracy of the ground 
state is thus given by : 



d g (X 



-^susy ) 



N/2 

E 

s=o 



N 
N/2 



(10) 



\ip )=Ae- r > s *\N/2,Q y ), 



(9) 



This collapse of the spectrum can be easily analyzed in 
the isotropic case 7 = 1 for which any state | S, M) is 
eigenstate with eigenvalue : 

E(S, M) = —j? (S(S + 1) - M 2 ) + A - 2hM. (11) 

For A < and h > 0, E(S, M) is minimum for M = S 
and at SUSY point E(S, S) = A for any S. 

The existence of this supersymmetric point enables us 
to locate the ground state for any value of the param- 
eters. Indeed, H describes a competition between the 
magnetic field h which aims to align the spin in the field 
direction, and the interaction terms which favors anti- 
ferromagnetic configurations. Thus, since at the SUSY 
point there exists a ground state lying in the maximum 
spin sector S = N/2, the ground state for h > h susy also 
lies in this sector and d g (h > h susy ) — D N / 2 — 1- Simi- 
larly, since singlet states |0,0) are also ground states at 
the SUSY point, they remain ground state for h < h susy 
and d g (h < h susy ) — Dq. In the thermodynamical limit 
(N — > 00), h susy goes to zero and the level crossing be- 
tween S — N/2 and S — ground states thus occurs at 
zero field. Nevertheless, for ft, = 0, it is clear that the 
ground states has a zero total spin for all 7 > 0. 

For odd N , the situation is more complex since the 
ground states of each spin sector do not degenerate as 
in the even N case except in the thermodynamical limit . 
In fact, as previously, the ground state belongs to the 
S = N/2 sector for h > h susy , and then switches to the 
other spin sector with decreasing S when h is lowered 
below h susy . Naturally, for h — 0, the ground state lies 
in the minimum spin sector and is given by all states 
1 1/2, 1/2). Its degeneracy D§ thus strongly depends on h. 
Nevertheless, in the thermodynamical limit, the region in 
which these level crossings occurs shrinks and converged 
to the zero field point so that the parity of N becomes 
irrelevant for the macroscopic physical quantities. 

To analyze the entanglement properties of the ground 
state, we focus on the concurrence C which has been 
introduced by Woofers to measure the two-spin en- 
tanglement. This quantity is obtained from the density 
matrix describing the state to be characterized. Here, 
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we concentrate on the ground state of H that can, as 
explained above, be degenerate. Thus, we must consider 
the thermal density matrix (at zero temperature) defined 
by: 



^ 0-9 



(12) 



"9 i= i 



where d g is the degeneracy of the ground state and where 
, d g } constitutes an orthogonal basis of the 



{\4>i),i = 1, 

dg-dimensional lowest energy subspace £q. Indeed, if we 
would consider the projector onto a specific state belong- 
ing to £q, the entanglement properties would strongly 
depend on this choice. Then, let p be the reduced den- 
sity matrix obtained by tracing out pth. over (N — 2) 
spins. Of course, in our system, the choice of the two 
spins kept is irrelevant because of the permutation sym- 
metry. Next, we introduce the spin-flipped density ma- 
trix p = a y ® o y p* o y ® Uy where p* is the complex 
conjugate of p. The concurrence C is then defined by: 



C = max {0, jjL\ — pi — p,3 — pa} ■ 



(13) 



where the pj are the square roots of the four real eigen- 
values of pp, classified in decreasing order. This con- 
currence vanishes for an unentangled two-body state 
whereas C = 1 for a maximally entangled one. As ex- 
plained in Ref. 0], it is further important to deal with 
a rescaled concurrence Cr = (N — 1)C to take into ac- 
count the coordination number of each spin. For a large 
number of spins, the difficulty comes from the trace step 
which requires operation in the full Hilbert space which 
is 2 w -dimensional. 

For h > h susy , this can be achieved since the ground 
state lies in the sector S = N/2. Indeed, (i) the 
S = N/2 subspace is nondegenerate so that the ther- 
mal density matrix of the unique ground state = 
J2 M a MW/2,M) simply writes p t h. = MW\ and (ii) 
the symmetry of the Dicke states \N/2,M) 24] allows 
us to write down in a simple form the reduced density 
matrix p in the standard basis {| |T), I tl)j I It); I II)}- 
One has 0: 

^ |2 (JV + 2Af)(JV + 2M-2) 



P22 = 



A I 



E 

M 



W(N- 1) 

,(N- 2M)(N + 2M) 
AN(N- 1) 



(15) 



^ l2 (N-2M)(N-2M-2) 
m = El^l W(N-l) (16) 

Pu = Yl a Ma*M+W(N + 2M)(N + 2M - 1) x 



y/(N - 2M + 1)(N- 2M + 2) 
4N(N- 1) ' 



(17) 



Furthermore, one has P23 = P33 = P22, and p* = p. 
The other matrix elements can also be computed for an 



arbitrary state belonging to S = N/2 subspace but since 
they vanish for the eigenstates of H because of (@J, we 
do not give them here. Note that for the Dicke states, 
one recovers the expressions given in Refs. [25lE^ . 

For h < h susy , the ground state lies in the minimum 
spin sector and is highly degenerate. The thermal density 
matrix is then simply the projector onto the subspace 
corresponding to S = for even N and S — 1/2 for 
odd N. The concurrence computed from such a density 
matrix is shown to vanish for N > 3 (2?| ■ 

Figure shows the behaviour of the rescaled concur- 
rence of the ground state as a function of the magnetic 
field. Note that for 7 = 1 and h > h susy the ground state 
is simply the Dicke state \N/2, N/2) for which C = 0. 
The same result also trivially holds in the large h limit 
for any 7. 




FIG. 2: Rescaled concurrence of the ground state as a function 
of the magnetic field for various anisotropy parameter 7 and 
for N — 10 3 spins (A = —1). Note that for any 7, one has 
Cr = at zero field. 



For h 



the worse situation is reached, at least 



for even N, since there is one ground state in each spin 
sector. As a consequence, the thermal density matrix is 
a sum of projector onto states of very different nature 
so that the trace operation is a rather difficult task for 
large N. Here, we have chosen to focus on \ipo) whose 
analytic expression Q allows one to compute its concur- 
rence for any even N. The coefficients olm entering in 
its decomposition onto the Dicke states are simply given 
by: 



a W/2-2j 



- v (N/2-2j) f N/2 \ 




(18) 



where we have set M = N/2 — 2j because of the symme- 
try Ili cr zl' ! /'o) — l^o)' Two limiting cases can be easily 
analyzed : the XY case (7 = 1) for which Cr = 0, and 
the Ising model (7 = 0) for which 77 = and Cr = 1. 
The rescaled concurrence of j^o) as a function of 7 is 
displayed in Fig. [3] 
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FIG. 3: Rescaled concurrence of the ground state \ipo) at the 
SUSY point as a function of the anisotropy parameter 7 for 
N = 10 4 spins. 

In the large N limit, the behaviour of the rescaled con- 
currence of l^o) can be computed and is given by : 

C R (l) = 1 - V7 (19) 

for < 7 < 1. Of course, this information does not 
enable us to conclude anything about the concurrence 
computed with the full thermal density matrix involving 
all spin sector but it certainly points out a nontrivial 
behaviour of the true ground state at h = h^ usy . 

Such a discontinuity of the concurrence at the transi- 
tion point has already been obtained in other frustrated 
spin models such as spin ladders or Heisenberg anti- 
ferromagnets in the Kagome lattice |l2| . However, in 
these systems the ground state entanglement properties 
are very simple (valence-bond like states) so that its 
concurrence is constant in each region of the phase dia- 
gram. If a jump of the ground state concurrence seems 
reasonable for a system displaying a first order quantum 
phase transition, it is not obvious that other measures of 
the entanglement would have shown discontinuity. More 
precisely, in the above mentioned example as well as in 
our model, the concurrence is found to become trivial 
for some parameter values. An interesting perspective 
would be to study other measures of the entanglement 
in systems displa ying a quantum phase transition such 
as the iV-tangle [28| or the Minkovskian-square norm 
of the Stokes tensor |2^| which investigate the N— spin 
entanglement. Indeed, the trace operation performed 
in the concurrence calculation undoubtedly kills some 
correlations between spins that could be captured by 
other types of measurement. 

We are indebted to B. Doucot and D. Mouhanna for 
fruitful and valuable discussions. 
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